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ABSTRACT 


Continuous Phase Frequency Shift Keying (CPFSK) is a 
potentially attractive modulation scheme with constant 
envelope and good spectral characteristics for energy 
constrained and band-limited satellite channels. 

This research deals with Orthogonal Frequency Shift 
Keying (FSK), Minimum Shift Keying (MSK) which is a special 
case of CPFSK, uncoded quaternary CPFSK and finally coded 
quaternary CPFSK. Orthogonal FSK is simulated by making the 
modulation index (h) equal to one, and all the other 
simulations are performed with h=1/2. 

A rate 1/2 convolutional encoder with constraint lengths 
(k), k=2, 3, and 4 are used in coded quaternary CPFSK 
simulations. Good coding gains are obtained with only a slight 
increase in receiver complexity. 

Soft decision with the Viterbi Algorithm was applied to 
all CPFSK and one MSK application. Hard decision was applied 


to Orthogonal FSK and another MSK application. 
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I. INTRODUCTION 


An increasing need for communication between individuals 
and nations is making satellite communication more and more 
attractive every day. Physical constraints of space based 
applications force us to choose bandwidth and energy 
efficient modulation techniques. 

Continuous Phase Modulation (CPM) is a digital modulation 
scheme with continuous phase and constant envelope. Continuous 
phase makes it spectrally efficient, and constant envelope 
seems desirable where non-linearities are present in space 
segment power amplifiers and repeaters. 

Coding is used to reduce reception errors. The trade-off 
in coding is a decrease in bandwidth efficiency. This thesis 
Mainly presents the research for rate 1/2 convolutional 
quaternary CPM schemes with different constraint lengths. 

The starting point for this research was the simulation of 
rate 2/3 trellis coded CPFSK. The complexity and apparent poor 
performance of this code forced us to first study orthogonal 
FSK, and then continue with MSK, uncoded quaternary CPFSK and 
finally finish with rate 1/2 trellis coded CPFSK. Due to the 
limitation of time, we were unable to return to the 2/3 rate 
codes. 


The simulations were done using PRO-MATLAB software run on 


Sun SPARC station computers. Hardware limitations and the 
complexity of the CPM schemes made this research difficult and 
very time consuming. 

Chapter II gives a brief review of Continuous Phase 
Modulation. Chapter III focuses on Orthogonal FSK, and two 
different receiver structures for MSK. Chapter IV presents 
uncoded quaternary CPFSK and coded quaternary CPFSK with 
different constraint lengths. Finally, Chapter V focuses on 
the final conclusions and opportunities for future research. 

An identical string of fifty thousand information bits has 
been used for each simulation to provide a comparison between 


systems. 


II. BACKGROUND 


This chapter focuses on the basic characteristics of 
Continuous Phase Frequency Shift Keying. In CPFSK, frequency 
shifts are used to transmit information. The phase of the 
CPFSK signal is continuous during frequency shifts. This 
prevents the large spectral side lobe effects that are 
observed in FSK, and PSK. There is also memory in CPFSK due to 
the continuous phase. Phase trellises will be presented in 


this chapter later on. 


A. CPFSK SIGNAL DESCRIPTION 


Equation (2.1) represents the transmitted signal in CPFSK, 


s(t) -|2Etcos 2ne.0+8 (6:1) +0, (Ze) 


where T,; is the symbol duration, E£, the energy per symbol, 


(Ref. 1] 


©(t;I) is the time varying phase modulation induced by the 
input data sequence, f, is the carrier frequency, and finally, 9, 
is the initial phase of the carrier assumed to be zero for 


simplicity. Equation (2.2) represents the phase of the carrier 


Pete interval nT. C< (mel)'T.. 


igi ab 
@(t;I)=nh) 1,+2nhgq(t-n7,) (2p 
k=-00 


The parameter h is called the modulation index. Typically his 


chosen to be a rational number oe where m and p are 


relatively prime positive numbers, in order that ©(t;I) takes 
on a finite number of states. The q(t) is a phase shaping 
pulse signal, I is the input signal. For M-ary CPFSK, M=2” 


and I, Us one of the (+1), 23) ,+(M-1)} information 


symbols. 


Phase shaping signals can be in a rectangular or raised 
cosine form. A rectangular form is used in this research. 
Phase trajectories for rectangular pulse shapes are piecewise 
linear. There are p phase states when mis even and 2p phase 
states when m is odd. Smoother phase trajectories can be 
obtained by using raised cosine pulse  shapes.Phase 


trajectories for binary CPFSK are illustrated in Fig. 2.1. 


B. ERROR PERFORMANCE 

Minimum Shift Keying (MSK) was used as a baseline 
modulation scheme; a special case of binary CPFSK with h=1/2. 
In coded CPFSK cases we also made comparisons with MSK, 
because, although we send two bits at atime, one of them is 
an information bit, and the other is a redundant coding bit. 


Equation (2.3) represents the probability of error in the 





Fig. 2.1: Phase trellis for binary CPFSK. 


presence of additive white Gaussian noise for MSK. [Ref. 2] 


p.=0(,|2—2) SS 


0 


The error performance can be improved by increasing the 
Signal energy, but in satellite applications there are 


limitations on power sources. 


Error performances for each modulation scheme are given at 


the end of each section. 


C. BANDWIDTH 

Bandwidth is the range of frequencies which contains 99 
percent of the total signal power. A smaller value of h will 
increase the bandwidth efficiency with an increase in 


Probability of exnuom 


III. FSK-MSK 
This chapter focuses on orthogonal FSK and MSK. The same 
modulator and bit mapper was used in all simulations. Input 


bits were mapped to modulation symbols according to Table 3.1. 


TABLE 3.1: MAPPING FOR ORTHOGONAL FSK AND MSK. 


inpuee bit mapper (1 
Che ~ CS alia =) al 
i = > 1 


A. ORTHOGONAL FSK 

Orthogonal FSK simulation is accomplished by using binary 
CPFSK with h=1. Coherent detection is used in the demodulator 
and hard decisions are made to detect the received signal. By 
making h=1, we have two phase states; 0 and a. The phase 
Paeeilromeer this system 15 illustrated in Fig 3.1. The odd 
symbols represent the system at state 0, and even symbols 
represent the system at state za. Multiplying the received 


pater ket 


Signal by e offsets the total phase to 0 at the end of 


each symbol and decreases the receiver complexity. The 
maximum value of the correlator outputs for the two modulation 


symbols is selected as the decision variable. 





Fig. 3.1: Phase trellis for Orthogonal FSK 


The MATLAB source code used for this system is listed in 
Appendix A. Equation (3.1) represents the probability of bit 


error for orthogonal FSK. [Ref. 2] 


ioe (3.1) 


The performance of this system is illustrated in Fig. 3.2. The 
orthogonal FSK Simulation was done fOr simulation 


verification, and as guidance to MSK simulation. 


5 Error Prabopiiity for Orthagonal FSK 
CoO —EeEe—e——e sss sss... SS a RCC oa Es oa oo 





Prob Of Error 


10-<¢ 
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Fig. 3.2: Performance of Orthogonal FSK 


B. MSK WITH VITERBI ALGORITHM DECODING 


MSK is a special case of binary CPFSK with h=1/2 where the 
phase changes +n/2 at every symbol (see Fig. 3.3). The total 
phase is +2/2 for odd symbols and the total phase is 0 or nm 
for even symbols. When the received signal is multiplied by 


eats it offsets the phase such that the total phase 


changes either by 0 or by zm at every symbol. This reduces the 


number of states required in the Viterbi decoder. [Ref. 3] 





Fig. 3.3: Phase Trellis for MSK 


The state trellis for the MSK waveform is illustrated in 
Fig. 3.4. There are two states, 0 and a. There are two paths 
out of each node corresponding to the two possible values of 
the input. Information bit 1 produces a change in state and 
information bit 0 produces no change. The Euclidean distances 
of the transmit symbol waveforms from the received signal are 
used as input to the Viterbi algorithm; the Viterbi algorithm 


chooses the path with the minimum total distance. 


10 


Fig. 3.4: State Trellis for MSK 


The transition matrix of the MSK with Viterbi algorithm 


decoding is illustrated in Table 3.2. 


TABLE 3.2: VITERBI TRANSITION MATRIX. 


AOR 21 83 
eee Ua 


Each row of the matrix corresponds to a state. The first 
column of the matrix represent the states that we are coming 
from, the second column represents the input value, and the 


third column represents the corresponding number of the point 


ial 


in the signal constellation. Columns 4, 5, 6 are the same 


respectively. The MATLAB source code for this system is given 


in Appendix B. 


Theory vs. simulation of the MSK bit error probabiliey 


ees 


with Viterbi detection is illustrated in Fig. 


Prob Of Error 





Fig. 3.5: Performance of MSK with Viterbi Detection. 


C. MSK WITH DETECTION USING TWO SYMBOLS AT A TIME 


Equation (2.1) can be rewritten as (3.2). [Ref. 4] 


a2 


| 
Suluter) “\ 2—tcos Ot wkcos Nagel 0). 
b (3E 2) 


2B sin[ @(t) pea (2% £ 2c) 


ie Se 


Equation (3.3) shows that the phase of an MSK signal increases 


or decreases linearly within the bit duration. [Ref. 4] 





@(t)=@(0)4n = ec) 


b 


The minus sign in (3.3) represents information signal 0, and 
the plus sign represents one. The signal phase can have values 


of +n/2 at odd multiple of T, or 0 and m at even multiple of 
eeecce @(oeemana Fig. 3.3). 


When (3.3) is inserted in (3.2), it is clear that the 


input sequence during (0,T,) affects only the quadrature 


component of (3.2), and not the in-phase component over the 


Paggealiang intemal (0;,27,). Similarly the input sequence 
@uring (72,22,) affects only the in-phase component of (3.2), 


and not the quadrature component over the signalling interval 


amet). 


Table 3.3 was constructed from the information discussed 


above. [Ref. 4] 


als 


TABLE 3.3: EFFECT OF INFORMATION BITS ON PHASE. 


INFORMATION BITS 





The hard decision in the demodulator is done according to 
Table 3.3, and the MATLAB source code for this system is 
illustrated in Appendix C. 

Theory vs. simulation of MSK bit error probability for 


detection using two symbols at a time is illustrated in Fig. 


326. 


Error Probopility for MSK 





Prob Of Error 


Fig. 3.6: Performance of MSK for detection using two 
symbols at a time. 


14 


IV. CODED CPFSK SCHEMES 


Channel coding can be used to reduce the number of errors 
when a message is transmitted over a noisy channel. To achieve 
this goal, redundancy is added to the information sequence. 

The encoder used in this research is a linear, sequential, 
finite state machine which is implemented by using shift 
registers. Rate 1/2 encoders with constraint lengths of 2, 3, 
and 4 are used. The transmitter block diagram for these 


systems is illustrated in Fig 4.1, and receiver block diagram 


mop illustrated in Fag. 4.2. 












CONVOLUTIONAL 
ENCODER 
RATE 1/2 













SOURCE MAPPER 


MODULATOR 


Fig. 4.1: Transmitter Block Diagram. 


Ss 















VITERBI 
DECODER 


SOURCE 
DECODER 





DEMODULATOR 


Fig. 4.2: Receiver block diagram. 


The bit mapping used in this chapter is illustrated in 


Table 4.1. Two coded bits are sent in each symbol duration. 


TABLE 4.1: MAPPING FOR CODED 4-CPFSK. 


input dibits mapper (I 
sa FB — 
Oly 23222 as = eco 
10 ---------- = ae 
11 <---<-<--=-<-- > 3 


The phase trellis for the uncoded quaternary CPFSK with 


16 


oe eee oe aerated in fig, 4.3. By combining coding 
techniques presented in this chapter, the probability of 
choosing the wrong information sequence will be reduced 


Significantly. 


9*pi/2 
8*pi/2 
7*pi/2 
6*pi/2 
5*pi/2 







Fig. 4.3: Phase trellis for quaternary CPFSK. 


The demodulator used in this chapter is illustrated in Fig 
4.4. As is seen, the received signal is first offset 
and then passed through four correlators. Equation (4.1) 
Tepeesciters thtecorkel atom ousput se 


ee 
Cae, ro (UW age ae Wot (4.1) 
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a 
Correlator 2 - 
= 


EUCLIDEAN 


DISTANCES 


ouad4dii>preo»peDeBF> vwvzeoo 


Correlator 4 


Fig. 4.4: Demodulator. 


Real parts of the correlator outputs are taken, and are used 
as elements of a four dimensional received signal vector. The 
euclidean distances of this vector are calculated in the 
comparator from each of the vectors in the signal 
constellation. The vectors of the signal constellation are 
determined as follows. 

In the absence of noise, the received signal will be 


r(t) =e 79" 2 (*k/Ts) Where the phase state @,=(0,m) and the 


transmitted information sequence J,=(-3,-1,1,3). Thus (gem. 


are eight vectors in the signal constellation determined by; 


18 


Cen= a Rel PeGees ec) dC] (4.2a) 
s 0 
4 Fe. es I.) 
eee fe ete Tete) (4.2b) 
+s e) 
= as peat. -2y) C/T, 
6 pei. a Rel f, e dt] (4.2c) 
Cumin (4.2d) 


Table 4.2 illustrates the Y,, in a tabulated form. Each row 


of Table 4.2 represents a single point in the signal 


constellation. 


TABLE 4.2: SIGNAL CONSTELLATION. 





Equation (4.1) is represented as (4.3) in the MATLAB 
program. 
TePOee 


GER SS ies ek Ve\ey (a in (43) 


k=0 


Le 


where At is the sampling period, and T,/At is the number of 


samples per symbol (N). In MATLAB simulation, Table 4.2 


becomes Table 4.3. 


TABLE 4.3: SIGNAL CONSTELLATION FOR MATLAB APPLICATION. 








A. CODED 4-CPFSK WITH CONSTRAINT LENGTH TWO 

The coding used in this section is created by using two 
shift-registers. [Ref. 5] The encoder is illustrated in Fig. 
4.5, and the state trellis for this system is illustrated in 
Figea oo 

The demodulator calculates the euclidean distances of the 
received signal vector from the eight vectors illustrated in 
Table 4.4. At each symbol duration, eight distance values are 
entered to the Viterbi algorithm and the true path is chosen 


according to minimum distance. 


20 


Most significant bit 


Input 


Least significant bit 


mig. 4.5: Encoder. 


9*pi/2 
8*pi/2 
7*pi/2 
6*pi/2 
5*pi/2 
4*pi/2 
3*pi/2 
2*pi/2 
1*pi/2 
0 
-1*pi/2 
-2*pi/2 
-3*pi/2 
-4*pi/2 
-5*pi/2 
-6*pi/2 
-T*pi/2 
-8*pi/2 
-9*pi/2 


States 





(a) 


(b) 


Fig. 4.6: Phase Trellis (a), and State Trellis for coded 4- 
CPFSK with k=2 (b). 


Za 


When comparing Fig. 4.3 and Fig 4.6 (a), it is clear Chee 
at each symbol duration, the probability of choosing the wrong 
State is reduced by this simple coding. Later, in the 
following sections, more complicated codes will be studied and 


the coding gains that are obtained will be seen more clearly. 


TABLE 4.4: SIGNAL CONSTELLATION. 


0,-1,0,—-N. GO, 1,0,N 1,020 eee eee 
O,=N,Cy=-l O;,N,0,1 <V,0,1,0 =Nat yee 








The transition matrix for this system is illustrated in Table 


oe he 


TABLE 4.5: VITERBI TRANSITION MATRIX. 


LO 1, 1 Ores ee ee 
2 5 1 2 6 Zee 2 oe 














Fach row of the matrix corresponds to a state. The first 
column of the matrix represent the states that we are coming 
from, the second column represents the input value, and the 
third column represents the corresponding point in the signal 
constellation that is illustrated in Table 4.4; the points are 
numbered rowwise. All other columns are the same respectively. 

The performance of this system is illustrated in Fig. 4.7. 
The coding gain obtained will not be worth what must be 


sacrificed in bandwidth efficiency (see Fig. 4.7). Since this 


PAP: 


was the initial experiment in the use of coding, it provided 
a baseline for our study. The MATLAB source code for this 


system is listed in Appendix D. 


4—CPFSK with ke2 


Error Probability for MSK ond 


Prob. OF Lrror 





Fig 4.7: Performance of the coded 4-CPFSK with k=2. 


B. CODED 4-CPFSK WITH CONSTRAINT LENGTH THREE. 

The coding used in this section was created by using three 
shift-registers. The encoder is illustrated in Fig. 4.8, and 
the state trellis for this system is illustrated in Fig 4.9. 
[Ref. 5] 

The demodulator calculates the euclidean distances of the 


received signal vector from the eight vectors presented in 


oe 


Most significant bit 





Least significant bit 


Fig 4.8: Encoder. 


States 
00 





Fig 4.9: State Trellis for coded 4-CPFSK with k=3. 
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Table 4.4. At each symbol duration, eight distance values are 
entered to the Viterbi decoder and the true path is chosen 
according to minimum distance. 

The transition matrix for this system is given in Table 
4.6. 


TABLE 4.6: VITERBI TRANSITION MATRIX. 


ee 2 822° 5 2.2 6 
Same os 6413 414 
Zim Ome 26 1.2 5 
Sy CU i. Se Se 8) ers ee oc a a 


The performance of this system is illustrated in Fig. 


ee O.. 


Error Probobility for MSK ond 4—CPFSK with k= 3 





Prob. Ot Error 
0 
’ 
> 


is) omelette! fii: PP Fri sass ses 
SS SRB SSH SS SBR DUBS SRG SSO OSS SSE SAD See 6 8 6 SSS a ae ee eee SC O05 Seinen eee eae ae ae moma 


Fig. 4.10: Performance of coded 4-CPFSK with k=3. 


Zo 


The coding starts working efficiently after 6.7 dB E,/N,: 


At this signal-to-noise ratio almost 1 dB gain is obtained and 


at 7.6 GB £,/N,, the coding gain reached almost 2 GB. 


Signal-to-noise ratios for 8.8 dB and 10 dB were also 
tested and no errors were detected in fifty thousand 
Informacion: bits. 

The MATLAB source code for this system is listed in 


Appendix E. 


C. CODED 4-CPFSK WITH CONSTRAINT LENGTH FOUR. 

The final coding studied, uses four shift registers. The 
encoder is illustrated in Fig. 4.11, and the state trellis for 
this system is illustrated in Fig. 4.12. [Ref. 5] 

In this application, the demodulator calculates the 
euclidean distances of the received signal vector from eight 
different vectors presented in Table 4.4. At each symbol 
duration eight distances are entered to the Viterbi decoder, 
and the true path is chosen according to minimum distance. 

The transition matrix of this system is illustrated in 
Table 4.7 and the MATLAB source code is listed in Appendix F. 

The performance of this coding is illustrated in Fig. 


4.13. The coding starts working efficiently at 5.9 dB £,/N,. 


At this signal-to-noise ratio almost 1 dB gain is obtained and 


at 6.8 dB £,/N,, the coding gain reached almost 2 GB. 
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Most significant bit 





Least significant bit 


Fig. 4.11: Encoder. 


States 





Oo N © oa HR W WN 


Fig. 4.12: Trellis for coded 4-CPFSK with k=4. 
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TABLE 4.7: VITERBI TRANSITION MATRIX. 


coe 


= 

PF OWN WN W © 
Pm DMNA MN ~I BR 
Jom Wwe ~~ TI MN Pf fe 
rP OWN WN W © 
WP TI HN oO HD Ww WNW 
aon FDNY OH WwW NWN 
Ww NM FP O FP OF N 
Jo mwbnhn £& NY WwW OD 
on FN DO HD W ND 
) [S) [as tS) eS) fF) 
@0—_.O}>.— {> __}- 3-3 


Signal-to-noise ratios of 7.7 dB, 8.8 dB, and 10 Geum. 
also tested and no errors were detected for fifty thousand 


infompmation bits 


Prob, Of Error 





Fig. 4.13: Performance of the coded 4-CPFSK with k=4. 
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V. CONCLUSIONS AND FURTHER RESEARCH 


This thesis mainly presented the simulation results for 
rate 1/2 convolutional quaternary CPM schemes with different 
constraint lengths. Using different codes, the minimum merge 
length increased, and larger minimum Euclidean distances were 
obtained. Better error performances were obtained as a result 
Sr these codes. 

This research can be extended in several ways. First, rate 
1/2 codes can be improved by using different constraint 
lengths, and shift register connections. The same encoders can 
also be tested with different modulation indexes. 

More complicated rate 1/2 codes with h=1/2, can be 
Simulated very easily by making small changes to the MATLAB 
source codes that are listed in Appendices D thru F. The 
Signal constellation that is presented in Table 4.4 can be 
used for every new code. 

Second, rate 2/3 codes for different constraint lengths, 
shift register connections, and modulation indexes can be 
studied, and best systems with minimum receiver complexity, 


and maximum power-bandwidth performances can be found. 


ENS) 


APPENDIX A 


This program creates am by n binary input matrix. 


Function [random Matrix] —sandeminy as 
a=rand (m,n) ; 

b=ones (m,n) *0.5; 
random_matrix=floor (a+b) ; 


This program is used as mapper. It converts 0 to -1 andl 


a of? 
ct 
O 
S 


function [mapper _output] =mapper (input) ; 
[m,a] =Size(input) ; 
for x=l:a 
y=ones (m,1)*2*(2-x) ; 
|r 2) ; 
end; 
mapper output=input.*z; 
Mapper output=mapper output’; 
mapper output=(mapper output) -1; 


This program is used as CPFSK modulator. Inputs to this 
program are "T=symbol duration", "h=modulation index", and 
"modulator _input=mapper output". 


a\P o\P ol? 


function [mod_output] =modulator(T,h,modulator_inpue); 
z=length (modulator input) ; 
tetad0=0; 
for sample=1:2z 
I=modulator_input (sample) ; 
1f sample== 
teta=0; *Sets the initial phase to zero. 
else 
teta0=teta0+modulator_ input (sample-1) ; 
teta=pi*h*tetad; 
end 
time=0:0.0005:T-0.0005; 
time=time+(sample-1)*T; 
xxX=pi*I*h* (time/T- (sample-1))+teta; 
mod_output (sample, :)=(exp(i*xx)); 
end 
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This program is used to add Gaussian noise to modulator 
GUGEUes 2NPUES PeuemIs programmeare "mM aS modulator output", 
"coeff as standard deviation of the noise", and "T as symbol 
duration". The constant 0.005 as sampling period. 


a? olP oP ol? 


function [awgn_noise] =awgn_ch(m,coeff,T) 

rand (*seed’ ,99); 

merida (’ normal’ ) ; 

awgn_noise=coeff* (rand(m,T/0.005) +j*rand(m,T/0.005) ); 


This program is used as CPFSK demodulator. Inputs to this 
program are; "T as symbol duration", "h as modulation index" 
"mod output as received signal". In this demodulator hard 
decision is used. 


a\? ol? ol? ol? 


function [demod_output] =demodulator(T,h,mod_output) ; 
Map=[(-1L lJw 
im, a] =saze(med soutput) ; 
@ame—-0:0.0005:T-0.0005; 
for sample=l1:m 
x=rem(sample, 2) ; 
phase_coeff=[(1-x)*pi]; 
phase=exp (i*phase coeff) ; 
Or Mmary=1:2 
xx=exp (i*time*pi*h*map(m_ary) /T) ; 
match=xx*conj (mod_output (sample,:)’); 
match=phase*matcnh; 
match_out=[match_out match] ; 
end 
[y, 1] =max(real (match_out)); 
demod_output (sample) =map(ceil(i/1)); 
match out=[ J; 
end 


cial 


functdon (Nb,1,¢e] =cheas 


% SYMBOL ERROR CHECK 

% Paul H. Moose 

% Naval Postgraduate School 
% 09-01-91 


+ This m-file locates the positions in vectors x and y that do 
noe 

¢ agree. It returns a one in e if they do not agree and a zero 
if they 

% do agree. e is the error vector if x and y are binary. I is 
a VeECreor 

% of error location numbers. Nb is the sum of the elements of 
e 

e 


=(x~=y) ; 
I=find(e); 
Nb=sum(e) ; 


% This program is used as master program. It calls different 
%* functions and simulates whole communication system. 


clear 

T=0.015; % "T" is used as symbol duration. 

h=1; * "h" is used as modulation index. 

m=50000- 

geal $ "m" and "n" are used to create binary matrix. 
coeff=0; % "coeff" is used as standard deviation of the noise. 


[input_matrix] =random(m,n) ; 
[mapper_output] =mapper (input matrix) ; 
[mod_output] =modulato(T,h,mapper_output) ; 
[awgn_noise] =awgn_ch(m,coeff,T); 
mod_output=mod_output+awgn_ noise; 
[demodulator output] =demodula(T,h,mod_output) ; 
error=check (mapper_output,demodulator_output) ; 


a2 


APPENDIX B 


This program is used as CPFSK demodulator. Inputs to this 
program are; "T as symbol duration", "h as modulation index" 
"mod output as received signal". In this demodulator hard 
decision is used. 


of olf ol ol? 


Punction [demod_ output] =demodulator(T,h,mod_output) ; 
map=[-1 1]; 
[m,a] =size(mod_output) ; 
fmme—0:0.0005:T; 
for sample=1:m 
mod=mod_output (sample, :)*exp(j*(sample-1)*pi/2); sOffsets 
the phase 
for m_ary=1:2 
xx=exp (i*time*pi*h*map(m_ary) /T) ; 
match=xx*conj (mod_output (sample,:)’); 
mMatch_out={match_out match] ; 
ena 
Gemod_ output (sample, :)=(real(match_out)); 
match_out=[ ]; 
end 
R= (j*demod_ output (:,1)+demod_output(:,2))/31; 
[demod_ output] =eucdis(2,R); ¢Calculates the euclidean distance 
Sof the received signal 


function D = eucdis(q,R) 
EUCLIDEAN DISTANCE METRICS 
Paul H. Moose 
Naval Postgraduate School 
06-17-92 


This M-file finds Euclidean distance of elements 

in vector R from 2*q unit amplitude vectors 

equally spaced on the unit circle. It stores these as rows 
GE 1. 


N=2°q; 

L=length(R) ; 

index=1:N; 

G@ph=2*pi/N; 

MO=exp (j*dph.* (index-1)); 

PO = 1: i 
D(1,:)=abs(R(1).*ones (MO) -MO) ; 
end 
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function PHN = softv (kK Ne, ae), 
Soft Viterbi Decoder 
Paul H. Moose 
Univ. degli Studi di Padova 
17-Os—en 


This M-file decodes k bit msgwords from 2%n real metrics 
(These may, for example, represent the "distance" of the 
received modulation value from each of 2*n modulation 
values.) 

* The state transition information for a 2°K state trellis ais 
alye 
% the 2*K by 3*2*k matrix T. Each of the 2°k entering paths 


IE, 
% each state has its source state (one of 2°K), path msgword 
( 
% of 2*k) and path codeword (one of 2*n) listed in the state 


row. 
S The path histories are kept in matrix PH that is 2°K by 
SND. 
% The path history for each state contains source state, path 
%* weight and path codeword for Np previous states. 
6 The output PHN is the update of PH, the new path history. 
% The decoded code word is in the last column of PHN. (They 
should 
6 "merge". 
% The past histories are undated on the basis of the 
"minimum 
% metric". You can change this to the "maximum metric" if 
desired as 
6 indicated in the comments in the code. 
Eons Ble 
K(j, 2) =D re 
Kai) saceae A 
X(j ] a) = 4 ’ 
rope oy 


wt = D(T(j,3*1))) +PH(T(373*1-2) 2 
if wt < X(j,2) %The < selects min metric 


)+PH(T(j,1),2); tpath weight 
path source state 
2) ;%path code word T(j,3) .Chg to T(j,2)aiaam 


° 
x) 


msgword. 


X(j,2)=wt; 
Kj, 1) = gee ene 
Ky, Seah Gola sChg to F(a, 208 for 
codeword. 
else 
end 
end 


* We now need to append old paths to new paths to get 
survivors. 

PHN(},:)= ([X(j,:) PH(X(j,1), 1233ipeeee 

end 


34 


% This program is used as master program. It calls different 
% functions and simulates whole communication system. All 
% other functions are same as in Appendix A. 


Clear 
Bemp={1 0 1 2.1.4;1 2 0.3);.%Is used as,Viterbi path 
matrix. 
f—O. 015; 
h=1/2; 
m=50000; 
n=1; 
eocer=0; 
[input_matrix] =random(m,n) ; 
[mapper output] =mapper (input_matrix) ; 
{mod_output] =modulato(T,h,mapper output) ; 
{awgn_noise] =awgn_ch(m,coeff,T); 
mod_output=mod_output+awgn_noise; 
[demodulator_output] =demodula(T,h,mod_output) ; 
viterbi path _matrix=temp; 
Tl=zeres (2,60) ; 
for x=l1:m 
D=demodulator output (x,:); 
eee eortv(1,1,20,TT,viterbdi path matrix,D) ; 
received signal (x) =TT(1,60) ; 
end 
Popout Matrix=input_matrix(19:m) ; 
received signal=received_signal (1:m-19) ; 
error=check (received_signal,input_matrix) ; 


f 


2)5) 


APPENDIX C 


This program is used as CPFSK demodulator. Inputs to this 


& 

fo] 

% program are; "T as symbol duration", "h as modulation index" 
% "in as received signal". In this demodulator hard decision 
° s 

6 1s used. 


function {demodulator output] =demodula(T,h, in) 

N=T/.0005; 

{k,a]=size(in); 

D =i:a, 

y3=-1; 

lene jee 
phil=exp (j*pi* (m-1) /2) *exp(j* (pi/2) .* (p/N)); 
phi2=exp (j*pi* (m-2) /2) *exp(j* (pi/2) .* (p/N) ) ; 
x=rem(m,2) ; 


cee (0) 
y=sum(real(in(@m-2,:)) . *abee( imag (57259. + 
sum(real(in(m-1,:)).*abs(imag(phil))); 
if ((y<0 & y3>0) | (y>0 &y3a<0)) 
a—as, 
else 
d=-1; 
end 
else 
y=sum(imag(in(m-2,:)).tabs(real( phi se + 


sum(imag(in(m-1,:)).*abs (real (phil) )); 
if ((y>0 & y3<0) | (y<0 & y3>0)) 


Vey 
demodulator output (m) =d; 
end 
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6 This program is used as master program. It calls different 
% functions and simulates whole communication system. All 
% other functions are same as in Appendix A. 


h=1/2; 

m=50000; 

n=1; 

coeff=0; 

[input matrix] =random(m,n) ; 

[mapper output] =mapper (input_matrix) ; 
{mod_output] =modulato(T,h,mapper_output) ; 
[awgn_noise]} =awgn_ch(m,coeff,T) ; 

mod output=mod_output+awgn noise; 
[demodulator_output] =demodula(T,h,mod_output) ; 
error=check (mapper output (2:m-2) ,demodulator_ output (4:m)); 


oy | 


APPENDIX D 


This program simulates the encoder that is given in 
Fig. 4.4. The input to this program is a me by 12 Dinas 
input matrix. At the end it gives a coded m by 2 binary 
output . First, two shift registers are initialized and then 
input are shifted. 


AP ol? o\® ol? ol? 


Function [tcm_output] =tem(anpuesnaeemes 

[m,a] =size (input matrix) ; 

ff1l=(O;input matrix(:,1)]; % Intializes the first £1ip eae 
ee Alpe aoe (alg gy elle 

fi2=— (Oy era); % Intializes the second flip-flop; 
fE2=fE2 (tee 

CemmoueDpie (:) 1) She, + Most significent bit. 

ECITGME DUE (=, 2)SEr2- * Least significicent bit. 


This program is used as mapper. It simulates the 
Table 4.1. 


al? ol? 


function {mapper output] =mapper (input) ; 
[m,a] =Size(input) ; 
fOr 2-134 
y=ones (m,1) *2* (3-x); 
z=(Z,yl; 
end 
Mapper output=input.*z; 
Mapper OQutPpUuE=mMapper CurEpUE F 
mapper_output=sum(mapper output) -3; 
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This program is used as CPFSK demodulator. Inputs to this 
program are; "T as symbol duration", "h as modulation index" 
"mod output as received signal". 


of al? ol? 


function [demod_output] =demodulator(T,h,mod_ output) ; 
Mapl=(-3 -1 3]; 
[m,a] =size(mod_output) ; 
pame=0>0.0005:T; 
for sample=l1:m 
mod=exp (j* (Sample-1) *pi/2)*mod_output (sample, :) ; 
for m_ary=1:4 
xx=exp (i*time*pi*h*mapl (m_ary) /T) ; 
mat chl=xx*coniMmod&) > 
Match_out=(match_out matchl']; 
end 
demod_ output (sample, :) =distance (real (match_out)); 
match out=[]; 
end 


This program calculates the euclidean distances of the 
received signal from signal constellation. 


ae oe 


function [dist] =distance(R) ; 


Mert rix=(O0 -1 0 -31 % This matrix represents the signal 
oe cal sconstellation. 
iOeod 0 
Si eh ae ik. 1G) 
esl) 0) - 1 
Ome O 1 
ale 01 0 
orl =O |S 
d=ones (8,1); 
G@a=a*R: 


dd=(dd-matrix) .*2; 
dist=sum(dd’ ); 


a9 


% This program is used as master program. It calls different 
6 functions and simulates whole communication system. All 
% other functions are same as in Appendix A. 


coeff=0; 
[input matrix] =random(m,n) ; 
[tem_output] =tcm(input matrisss 
[mapper output] =mapper (tem outpue); 
{mod_output] =modulato(T,h,map) ; 
[awgn_ noise] =awgn_ch(m,coeff,T) ; 
mod_output=mod_output+awgn noise; 
[demodulator output] =«demodula(T, i, mod omepue :, 
viterbi_path_matrix=temp; 
THhazeros (27/220): 
£Onex= ikem 
D=demodulatorseutpul ea); 
[TT] =softv(2,1,40,TIT,viterbi path mame, b); 
received signal (x) =TT(1,120) ; 
end 
ww=mb (2, received signal) ; 
WW=ww(1:2:2*m) ; 
error=check (input_matrix(1:m-41)/’,ww(42:m)); 
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APPENDIX E 


This program simulates the encoder that is given in 
PaG-ea.. ihe input COuEhniS program”is a m by”1 binary 
input matrix. At the end it gives a coded m by 2 binary 
output . First, the three shift registers are initialized 
and then input is shifted. 


AP olf olf olf ol? 


mum@eGELOn [tcm output] =tem(input_ matrix) ; 
{m,a] =Ssize(input matrix) ; 


meee O; input maltrix(e)1)); * Initialize the first flip-flop. 
motel | lesm, 1) ; 
fro? = (10; fF) ; % Initialize the second flip-flop. 
m= feo Wem, 1) 
fete [Oot £2); % Initialize the third flip-flop. 


ees & 3)( 1 > Mizu) ; 
tem output (:,1) =abs (abs (f£1-f££2) -££3) ; + Most significant bit. 
Bemmeout pute 2) =££2; % Least significant bit. 
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function [PHN]= viterbi(k) Ne, ee ies, 
% Soft Viterbi Decoder 
% Paul H. Moose 
% Mercury Digital Communications 
% 06-09-91 
% This M-file decodes k bit msgwords from 2%n real metrics 
% (These may, for example, represent the "distance" of the 
% received modulation value from each of 2*n modulation 
values.) 
% The state transition information for a 2*K state trellis is 
in 

the 2*K by 3*2*k matrix T. Each of the 2*k entering paths 


|e 
% each state has its source state (one of 2*K), path msgword 
( 
% of 2*k) and path codeword (one of 2%n) listed in the state 


% The path histories are kept in matrix PH that is 2°K by 
3*Np. 

% The path history for each state contains source state, path 
%* weight and path codeword for Np previous states. 

% The output PHN is the update of PH, the new path history. 
%* The decoded msg word is in the last column of PHN. (They 
should 

6 "merge". 

% The past histories are undated on the basis of the 
"minimum 

6 metric". You can change this to the "maximum metric" if 


P=Pn se 

wt=T(:,3:3:3*2°k)'; 

UX=T (5 1:323*2"> k-=2) "+ 

aa=D(wt(:)); 

bbD= Bee ee) a 

wt(:)=saa+bb; %tThis contains all weights(columns) for each 
state (row) 

[a,b]=min(wt); * Use max(wt) here for maximum 

X(:j,2) 2a"; 


X (2, 1) =diag CT (253-562 

X(:,3)=diag(T(:,3.*b-1));spath msgword. Chg to 3.*b jem 
codeword 

sWe now need to append old paths to new paths to get 
Survivors. 

PHN=(X PH(X(:,1),1:3*Noaaiel 
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6 This program is used as master program. It calls different 
% functions and simulates whole communication system. All 
% other functions are same as in Appendix A. 


iW 
Pow Oo 
DN WP 
WN WH 
FOWO 
WH ON 
Db HOB tb 
WN N 
Own 


T=0.015; 
h=1/2; 
m=50000; 
aol. 
coeff=0; 
[input_matrix] =random(m,n) ; 
[tcm_output] =tcm(input_matrix) ; 
[mapper_output] =mapper (tcm_output) ; 
[mod_output] =modulato(T,h,map) ; 
[awgn_noise] =awgn_ch(m,coeff,T) ; 
mod _output=mod_output+awgn_noise; 
[demodulator_output] =demodula(T,h,mod_output) ; 
viterbi_path_matrix=temp; 
TT=zeros (4,120) ; 
for x=1:m 
D=demodulator_output (x,:); 
eevee roa(2,240,TT,viterbi path matrix,D); 
received signal (x) =TT(1,120) ; 
end 
ww=mb (2,received_signal) ; 
ww=ww(1:2:2*m) ; 
error=check (input (1:m-41)’,ww(42:m) ); 
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APPENDIX F 


This program simulates the encoder that is given in 
Fig. 4. The input to this program@is aem Dy sera y 
input matrix. At the end it gives a coded m by 2 binary 
output . First, the four shift registers are initialized 
and then input is shifted. 


AP ol olP olP o\P 


function [tcm_ output ] =ceminpurenaeais = 
[m,a] =size (input_matrix) ; 

fil= (O07 tnpucumacr we(e ee, 

EG dea eee si, ele 

£E2= (0; Eee 

Eel See Gee m,. 1): 

£i3=[(0-f£E2) - 

ES ee See: Meee 

EE4= (07 taal; 
ff4=ff4(1:m,1) 
ECM Out pure 


,1)=abs (ff1-abs (£FF2=£f4) ) ; 
ECMACUEBIIE:,; 2) — 


EES; 


4.4. 


This program is used as master program. It calls different 
functions and simulates whole communication system. All 
other functions are same as in Appendix A. 


oP ol? ol? 


clear 
temp= [ 


ab 
4 
5 
7 
al 
3 
5 


POWNWNW © 
BNO UD UTR 
Nowra U BPR 
rFPOWNWNW © 
Wr ~V0O OO @MON 
OnPN WOW ND 
WNRrROrFORN 
“J O1rW DN BN W OD 
OnP NY ®W AWN 
WI iS) (SO) (SS le bk 
OnPFrF UF BUM 


seo O 15; 
h=1/2; 
m=50000; 
n=1; 
coeff=0; 
[input_matrix] =random(m,n) ; 
[tem_output] =tcm(input_matrix) ; 
[mapper_output] =mapper (tcm_output) ; 
[mod_output] =modulato(T,h,map) ; 
[awgn_noise] =awgn_ch(m,coeff,T); 
mod_output=mod_output+awgn noise; 
[demodulator output] =demodula(T,h,mod_output) ; 
viterbi_path_matrix=temp; 
TT=zeros (8,120) ; 
for x=1:m 
D=demodulator output (x,:); 
[TT] =viterbi(2,40,TT,viterbi_ path _matrix,D); 
received signal (x) =TT(1,120); 
end 
ww=mb (2, received signal) ; 
ww=ww(1:2:m*2); 
error=check (input_matrix(1:m-42)',ww(43:m)); 
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